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2.  Producing  Quantitative  Time-Series  of 
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Accomplishments/New  Findings: 

1  Original  Goals 

The  stated  goals  for  this  project  can  be  summarized  as  the  following: 

1.  Develop  a  method  to  inter-compare  the  flare-predictive  capability  of  parameters  derived  from 
photospheric  vector  magnetic  field  data. 

2.  Apply  the  Minimum  Coronal  Current”  algorithm  [23]  to  photospheric  vector  magnetic  field 
data  with  the  goal  of  quantifying  the  coronal  magnetic  complexity  as  it  relates  to  flare  produc¬ 
tivity,  and 

3.  Acquire  and  prepare  vector  magnetogram  data  to  test  and  demonstrate  the  feasibility  of  the 
above  methods,  ensuring  that  the  data  are  sufficiently  quantitative  and  accurate  as  to  be  ap¬ 
propriate  for  this  project. 

In  short,  we  have  full-filled  our  commitment;  as  described  below,  in  some  cases  the  goals  were 
hindered  by  unforeseen  problems,  yet  in  others  our  initial  goal  has  been  surpassed.  In  both  situa¬ 
tions,  we  continue  our  close  collaboration  with  scientists  at  many  institutions,  including  the  National 
Center  for  Atmospheric  Research  (Drs.  S.  Gibson  and  Y.  Fan),  U.  C.  Berkeley  (Dr.  G.  Fisher/the 
MURI  team),  Montana  State  U.  (Drs.  D.  Longcope,  R.C.  Canfield  and  other  MURI  team  members). 
Below  we  summarize  our  accomplishments  and  new  findings  from  this  project. 

2  Producing  Quantitative  Time-Series  of 
Photospheric  Vector  Magnetic  Field 

A  significant  amount  of  the  effort  for  this  project  has  been  directed  to  producing  time-series  of 
photospheric  vector  field  data  that  are  adequately  characterized.  Both  and  random  uncertainties 
must  be  understood,  quantified,  and  included  such  that  variations  observed  in  the  time-sequence 
are  trusted  to  educate  us  on  the  solar  condition. 

2.1  The  General  Data  Reduction 

The  data  used  in  this  project  are  from  the  U.  Hawai‘i/Mees  Solar  Observatory  Imaging  Vector  Mag¬ 
netograph  at  Haleakala  [28,  13,  12].  Briefly,  the  spar-mounted  IVM  was  designed  to  minimize  both 
instrumental  polarization  and  internal  turbulence.  A  four-frame  polarization-modulation  sequence 
is  employed,  using  a  Fabry  Perot  etalon  to  sample  30  wavelength  positions  across  the  magnetically 
sensitive  (geff  =  2.5)  Fel  630.25nm  spectral  line  in  under  2  minutes. 

Corrections  to  the  raw  data  are  performed  to  remove  spatial  and  polarization  distortions  from 
both  the  telescope  system  and  atmospheric  seeing  [28,  13].  The  data  are  then  demodulated  to 
produce  Stokes  spectra  ([/,  Q,  U,  V])  at  each  pixel  [13].  As  a  final  step,  the  spectra  are  binned  to 
2562  1.1"  pixels.  The  final  polarization  noise  (normalized  by  I_continuum)  is  of  order  2  x  10~3  in 
these  data. 

To  derive  the  magnetic  flux  vector  from  the  resulting  spectropolarimetric  data,  we  employ  a 
forward-integration  scheme  based  upon  the  equations  of  [14, 15],  which  produces  results  comparable 
to  a  full  least-squares  inversion  [12],  The  resulting  maps  of  the  observed  line-of-sight  and  transverse 
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magnetic  components  have  uncertainties  of  order  10G  (Bg),  25G  (Bt),  with  an  uncertainty  in  the 
azimuthal  angle  of  approximately  5  deg. 

The  180deg  ambiguity  in  the  observed  transverse  component  is  resolved  using  an  automated 
iterative  procedure  which  minimizes  the  difference  between  the  observed  field  and  a  force-free  field 
computed  using  the  Bz  component  and  which  additionally  minimizes  the  field’s  divergence  [4]. 
With  a  spatial  map  of  the  heliographic  Bx,By,Bz  components  we  derive  corresponding  maps  of 
their  uncertainties,  computed  using  the  noise  in  the  observed  fields,  the  uncertainties  returned  from 
the  inversion,  and  then  the  coordinate  transforms  and  error  propagation  as  described  in  [21]. 

The  temporal  sequence  of  vector  magnetic  field  maps  are  then  aligned  on  a  sub-pixel  grid  for 
pointing  variations  and  trimmed  for  edge-data  artifacts.  The  final  data-cube  consists  of  [a:,  y,  B,  t] 
and  is  stored  in  an  easily  accessible  “structure”  format  of  the  Interactive  Data  Language  system 
for  analysis. 

2.2  Recently  Implemented  Improvements  and  Enhancements 

Over  the  course  of  this  contract,  the  P.I.  designed  and  implemented  improvements  in  the  following 
areas: 

1.  Flat  Field  Corrections:  The  flat-fielding  algorithm  was  updated  to  account  for  the  temporal 
variability  of  the  observed  fringe  patterns;  it  also  now  determines  the  excess  scattered  light. 

2.  De-Blur/De-Stretch:  This  algorithm  was  improved  both  for  speed  and  to  avoid  some  spurious 
noise  issues. 

3.  Data-Handling:  Algorithms  and  procedures  were  developed  (see  appendix)  to  more  effi¬ 
ciently  and  accurately  co-align,  trim,  store,  and  analyze  large  numbers  of  magnetograms 
in  a  platform-independent  manner. 

2.3  Algorithm  to  Account  for  Seeing  Effects: 

Leka  &  Rangarajan,  Solar  Phys .,  2001. 

The  distortions  in  ground-based  data  caused  by  turbulence  in  the  Earth’s  atmosphere  are,  to  a 
large  extent,  minimized  by  the  de-blur/de-stretch  algorithms  included  in  the  IVM  data  reduction 
package.  The  overall  quality  of  the  atmospheric  seeing  can  significantly  change  the  inferred  magnetic 
vector  and  all  quantities  derived  from  it  (Figure  1).  A  quantitative  study  of  the  magnitude  and 
nature  of  these  changes  was  undertaken  [20].  The  primary  results  included  a  surprising  degree  of 
variation  in  derived  quantities  (e.g.,  the  magnitude  of  the  vertical  electric  current  Jz,  the  degree  of 
magnetic  shear,  etc.)  due  to  modest  variation  in  seeing,  and  additionally  that  the  variations  can 
be  region-specific  (i.e.,  the  results  are  not  generalizable  between  data  from  different  active  regions). 
A  method  is  described  for  first  determining,  and  then  accounting  for,  the  seeing-induced  variations 
in  the  uncertainties  for  the  raw  data  and  all  quantities  derived  from  them. 
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1.0  0.6  0.2 
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Figure  1:  An  example  of  the  effects  of  worsening  seeing  on  one  typical  parameters  derived  from  vector 

magnetic  field  data.  The  total  signed  magnetic  flux  (left)  for  two  active  regions,  AR7981  ( - )  and  AR8771 

^ - ))  positive-polarity  (top)  and  negative-polarity  (bottom)  axe  plotted  as  a  function  of  a  relative  seeing 

measure  <Sr,  which  is  simply  the  ratio  of  the  blurred  data  to  the  reference  data  (i.e.,  *Sr=1  is  the  original 
dataset,  «Sr  <  1  is  worsening  seeing.  In  the  case  of  AR8771  <Sr=0.2  resulted  when  a  4.4"  blur  was  applied  to 
the  raw  images  while  for  AR7981  the  same  applied  blur  resulted  in  a  slightly  worse  <SR.)  The  error  bars  for 
the  data  points  are  shown,  as  is  a  linear  fit  to  the  effects  of  worsening  seeing.  For  direct  comparison  between 
the  two  active  regions,  the  effects  are  plotted  as  a  fractional  change  (right).  Here,  the  fit  is  plotted,  and  the 
error  bars  reflect  the  standard  deviations  to  the  fit.  (From  [20].) 


3  Parameterization  of  Active  Region  Magnetic  Distribution: 

Leka  &  Barnes  Astrophys.  J.  2003a. 

3.1  Abstract 

Using  photospheric  vector  magnetic  field  data  from  the  U.  Hawai'i  Imaging  Vector  Magnetograph 
with  good  spatial  and  temporal  sampling,  we  study  the  question  of  identifying  signatures  unique 
to  flare  events  in  parameters  derived  from  B.  In  this  first  of  a  series  of  papers  we  present  the 
data  analysis  procedure,  and  present  sample  results  focusing  only  on  three  active  regions  (NOAA 
Active  Regions  #8636,  #8771,  and  #0030),  two  of  which  produced  three  flares  (two  M-class  and 
one  X-class),  as  compared  to  (most  importantly)  a  flare-quiet  epoch  in  a  third  comparable  region. 
Quantities  such  as  the  distribution  of  the  field,  the  distribution  of  the  spatial  gradients  of  the  field, 
vertical  current,  current  helicity,  “twist”  parameter  a  and  magnetic  shear  angles,  are  parameterized 
and  examined  for  differences  in  overall  magnitude  and  evolution  between  the  flare  and  flare-quiet 
examples.  The  variations  expected  due  to  atmospheric  seeing  changes  are  explicitly  included.  In 
this  qualitative,  “single  quantity  considered  at  a  time”  approach  we  find  (1)  no  obvious  flare- 
imminent  signatures  from  the  plain  magnetic  flux  vector,  the  total  magnetic  flux,  the  total  current 
current,  or  free  energy  quantities;  (2)  we  find  counter-intuitive  but  distinct  flare-quiet  implications 
from  the  mean  spatial  gradient  of  the  field,  the  chirality-term  of  the  vertical  current,  and  the 
kurtosis  of  the  distribution  of  magnetic  shear  angles;  (3)  we  find  flare-specific  or  flare-productivity 
signatures,  although  sometimes  weak,  using  the  kurtosis  of  the  fields’  spatial  gradients,  the  mean 
vertical  current  and  current-helicity  densities,  mean  magnetic  shear  angles,  and  net  current.  The 
strongest  results  are,  however,  that  (4)  in  ensuring  a  flare-unique  signature,  numerous  potential 
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Figure  2:  Images  of  continuum  (top)  and  Bz  (bottom,  with  ±  100  G  contours)  of  NOAA  AR0030  (left), 
AR8636  (middle)  and  AR8891  (right).  Axes  are  approximately  in  Mm,  and  black  triangles  are  masked-out 
field  stops  (adapted  from  [18]) 


candidate  signatures  (in  both  the  variation  and  the  overall  magnitude  of  a  parameter)  are  nullified 
on  account  of  their  behavior  in  a  flare-quiet  region,  and  hence  (5)  this  qualitative  “single  quantity 
considered  at  a  time”  approach  is  inadequate.  To  address  this  latter  point  we  employ  a  quantitative 
statistical  approach  in  Paper  II  [19]. 

3.2  Example:  Magnetic  Flux  and  its  distribution 

The  total  magnetic  flux  is  a  quantitative  measure  of  an  active  region’s  size  which  is  well  correlated 
with  its  productivity  for  energetic  events  [7,  27,  5,  31].  Using  the  magnetic  flux  as  a  measure  of  size 
rather  than  e.g.,  the  white-light  area  provides  a  more  physical  description  of  the  energy  available 
for  such  events.  We  consider  here  the  total  unsigned  flux  $tot  =  X)  |-BZ|  dA  as  well  as  the  unsigned 
total  flux  for  each  polarity  separately  4>+  =  JX-B*  >  0)  dA,  =  Y1(BZ  <  0)  dA. 

We  also  consider  the  net  flux  imbalance  |4>„et|  =  |  dA\,  a  quantity  which  has  been  asso¬ 
ciated  with  flare  activity  [34,  31],  although  with  little  discussion  as  to  the  physical  influence  any 
(local)  imbalance  may  have.  Nonetheless  we  consider  this  quantity  but  take  its  absolute  value  to 
avoid  hemispheric  biases  imparted  by  the  asymmetric  nature  of  the  spatial  flux  distribution  between 
preceding  and  following  polarities. 

Sample  results  for  three  active  regions  (of  the  full  data  considered,  see  Table  2), for  total  and 
net  flux  are  shown  in  Figure  3.  The  three  regions  considered  are  comparable  in  size  (4-6  xlO22 
Mx),  although  flare-quiet  AR8891  is  the  largest.  None  of  the  regions  is  flux-balanced,  AR8636  the 
farthest  from  balance  by  1022  Mx;  this  is  hardly  a  surprising  result  given  a  restricted  (even  if  large) 
field-of-view. 

Both  flaring  regions  show  small  variations  on  the  timescales  shown  here,  with  some  flux  growth 
in  AR8636.  However,  the  changes  are  generally  not  beyond  the  uncertainties,  especially  when  seeing 
is  accounted  for. 

The  distribution  of  the  components  of  the  magnetic  field  (Bz  and  Bh)  over  the  active  region  is 
quantitatively  described  by  the  statistical  quantities  mean,  standard  deviation,  skew,  and  kurtosis. 
We  consider  the  distributions  of  the  spatially  sampled  (denoted  by  “s”  rather  than  “(x,  y)”  for 
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simplicity)  Bz(s),  Bh(s)  =  ]/Bx(s)2  +  By(s)2,  and  B{s)  =  ^Bx{s)2  +  By{s)2  +  Bz(s)2,  the  latter 
two  being  positive-definite  quantities. 

We  present  here  the  mean  of  the  B(s)  distribution,  the  skew  of  the  Bz(s)  distribution,  and 
the  skew  of  the  Bh(s)  distribution  (Fig.  3).  The  mean  pixel-averaged  magnetic  flux  similar  for  all 
three  regions  (400G-600G);  while  there  is  variation  in  this  quantity  prior  to  all  flares  shown  here, 
there  is  similar  variation  during  the  observations  of  flare-quiet  AR8891.  The  ability  of  skewness 
discriminate  otherwise  cloaked  aspects  of  the  fields’  distribution’s  is  clear  when  comparing  AR0030 
and  AR8891:  the  skew  of  Bz(s)  is  positive  for  the  former  but  negative  for  the  latter  allowing  us  to 
identify  a  significant  contribution  from  the  negative  tail  in  the  latter;  this  is  consistent  with  the  fact 
that  the  (signed)  net  magnetic  flux  (not  shown)  for  AR8891  is,  in  fact,  negative,  skewness  shows 
some  variation  in  all  three  regions  (a  marked  variation  in  AR8636  although  with  similar  changes 
due  to  the  seeing  variations).  The  overall  skewness  level  of  Bh(s)  is  similar  for  all  three  regions: 
it  is  positive  indicating  a  contribution  from  a  high-field-strength  tail,  i.e.,  penumbral  fields.  The 
skewness  increases  prior  to  the  M-flare  in  AR0030,  shows  an  overall  decrease  before  the  flare  in 
AR8636,  and  generally  increases  in  AR8891. 

In  summary,  using  the  distributions  of  the  magnetic  flux  vector  for  these  examples,  we  find 
no  obvious  and  consistent  difference  between  the  flare-productive  and  flare-quiet  regions,  and  no 
obvious  and  consistent  evolution  that  occurs  in  the  pre-flare  periods. 

3.3  Example:  The  Magnetic  Helicity  Density  hc 

To  describe  the  non-potential  nature  of  the  active  region  magnetic  fields  one  may  use  the  magnetic 
helicity  [2];  a  conserved  quantity,  the  describes  the  linkages  of  magnetic  systems.  In  practice,  a 
number  of  caveats  (including  the  possible  forced  state  of  the  solar  photosphere  and  the  limitation 
of  these  data  to  sampling  a  single  height  in  the  solar  atmosphere)  limits  what  we  can  immediately 
derive  to  the  vertical  component  of  the  current  helicity  density: 

hc(s )  =  Bz(s)VhxBh  (1) 

=  (2> 

This  quantity  is  determined  spatially  and  as  such  we  can  examine  the  moments  of  its  distribu¬ 
tion.  Additionally  we  compute  the  total  (unsigned)  He, tot  —  ]C  l^c  dA\  and  net  (signed)  \HC)flet\  = 

|  X)  hcdA\  current  helicity  over  the  active  regions,  the  latter  being  akin  to  the  current  helicity  im¬ 
balance  described  in  [2].  In  Figure  4  the  temporal  variation  of  these  quantities  as  well  as  the  mean 
of  the  current  helicity  density  distribution  hc. 

Interestingly,  we  find  that  overall  the  largest  total  current  helicity  occurs  in  AR8891  which  is 
the  flare-quiet  region  during  this  observing  epoch.  All  three  regions  show  temporal  variations  and 
one  could  argue  there  is  a  decrease  in  the  total  current  helicity  prior  to  the  flares;  there  is  a  similar 
decrease  during  the  observations  in  AR8891,  indicating  that  this  variation  is  not  unique  to  the 
flaring  condition.  On  the  other  hand,  the  net  (signed)  current  helicity  is  significantly  different  from 
zero  for  both  AR0030  and  AR8636,  while  it  is  consistent  with  zero  for  the  entire  observing  period 
for  AR8891.  This  is  consistent  with  some  examples  in  [2]  for  their  “imbalance”  parameter. 

Also  consistent  with  [2]  we  find  a  larger  variability  in  hc  for  the  flaring  regions  than  with  the 
flare-quiet  AR8891,  as  well  as  a  larger  overall  level  of  this  quantity.  Additionally,  for  both  M- 
class  flares  one  could  argue  a  decrease  in  hc  occurs  prior  to  the  flare  (although  it  is  not  clear  such  a 
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Figure  3:  Parameters  discussed  in  the  text  for  (left: right)  AR0030,  AR8636,  and  AR8891.  The  start-times 
of  the  flares  as  determined  by  the  GOES  SXR  light  curve  are  indicated  by  vertical  grey  lines,  the  X3.0  (thick) 
and  M1.8  flares  in  ARD030,  and  the  M1.1ERU  flare  in  AR8636;  AR8891  did  not  produce  any  flares  during 
this  epoch.  The  x-axes  indicate  the  UT  time,  y-axes  are  in  the  relevant  units,  and  1-a  error  bars  are  plotted 
(for  the  3-<r  data),  including  the  expected  variation  due  to  measured  changes  in  the  seeing  conditions  (the 
asymmetric  error  bars).  Shown  (top:bottom)  are  (a)  Total  unsigned  signed  magnetic  flux,  (b)  the  absolute 
value  of  the  net  magnetic  flux  imbalance,  (c)  The  mean  (pixel-averaged)  magnetic  flux  magnitude,  (d)  the 
kurtosis  of  the  Bz(s)  distribution,  and  (e)  the  skew  of  the  Bh(s)  distribution  (adapted  from  [18]). 


decrease  occurs  prior  to  the  X-class  flare)  where  no  decrease  of  similar  magnitude  occurs  in  AR8891. 
A  similar  decrease  in  the  magnitude  of  hc  is  observed  by  [22]  prior  to  an  X-class  flare  in  localized 
portions  of  a  different  active  region. 

To  summarize,  the  magnetic  helicity  displays  a  number  of  event-unique  signatures,  including  a 
small  HcM,  significant  .He, net  and  7£,  and  a  clear  decrease  in  hc  prior  to  two  of  the  three  flares. 

3.4  Summary 

This  and  our  subsequent  investigations  begin  with  the  null  hypothesis,  that  there  is  no  detectable 
signature  of  an  impending  energetic  event.  As  such  we  have  selected  data  to  specifically  not  be 
biased  for  flaring-only  regions  or  epochs,  but  to  include  (1)  flaring  epochs  and  (2)  and  flare-quiet 
epochs  from  the  same  region,  even  on  the  same  day,  and  (2)  flare-quiet  regions  which  had  been 
given  a  high  probability  of  flaring  (cf.  Table  2).  Without  attempting  to  indict  the  null  hypothesis 
it  is  impossible  to  determine  whether  there  is  a  unique  situation  in  the  solar  atmosphere  which 
produces  energetic  events.  In  other  words,  for  signature(s)  to  be  uniquely  related  to  solar  energetic 
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Figure  4:  Same  as  Fig.  3,  but  for  (a)  the  total  (unsigned)  current  helidty,  (b)  the  net  (signed)  current 
helicity,  (c)  the  mean,  (d)  skew,  and  (e)  kurtosis  of  the  current  helidty  density  distribution. 


events,  they  must  also  not  be  present  at  times  when  no  energetic  event  is  produced. 

Three  additional  issues  must  be  reiterated  here.  First,  we  explicitly  indicate  how  changes  in 
terrestrial  atmospheric  conditions  can  influence  the  results  (e.g.,  Figs.  3,  4).  One  simply  cannot 
make  ad  hoc  assumptions  about  the  variations  of  the  measured  quantities  using  data  which  were 
themselves  obtained  under  different  conditions  (both  solar  and  terrestrial). 

Second,  all  quantities  considered  here  (cf.  Table  1)  are  based  on  the  physical  situation  at  the 
photosphere.  Observational  biases  have  been  removed  (e.g.,  the  use  of  heliographio-plane  rather 
than  observational-plane  magnetic  quantities,  and  the  use  of  unsigned  quantities  where  appropriate 
to  avoid  hemispheric  biases),  or  acknowledged  (the  effect  of  any  limitation  on  the  field  of  view  of 
the  instrument,  no  matter  now  large,  on  “net”  quantities,  for  example). 

Third,  by  sampling  the  active  regions  both  temporally  and  spatially,  and  then  using  the  mo¬ 
ments  (and  where  appropriate,  the  total/net  quantities)  of  the  resulting  spatial  distributions,  we 
parameterize  the  magnetic  state  of  the  photosphere  and  quantitatively  allow  for  uncertainties  in 
the  data  and  observing  conditions.  In  this  manner,  we  aim  to  detect  subtle  variations  caused  by 
the  evolving  state  of  the  photospheric  magnetic  field  and  avoid  subjective  “by-hand”  examination 
of  localized  changes  which  might  occur;  the  latter  approach  is  crucial  when  (as  we  have  now),  the 
data  approaches  a  large,  statistically  significant  quantity. 

In  summary  we  find  no  obvious  flare-event  signatures  using  the  parameterization  of  the  magnetic 
flux  vector  distribution,  the  total  current,  or  free  energy  quantities. 

Counter-intuitively,  we  find  that  a  broad  distribution  of  inclination  angles  is  unique  to  a  flare- 
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quiet  episode,  as  is  a  greater  mean  magnitude  of  the  spatial  gradient  of  field-strength.  The  chirality- 
related  component  of  the  vertical  current  I^t  is  larger  overall  for  the  flare-quiet  region,  as  is  the 
kurtosis  of  the  shear-angle  distributions. 

Flare-specific  signatures  include  weak  indications  of  increased  variability  and  possibly  a  pre¬ 
event  rise  in  KflVfc-BI),  increased  variability  in  J®et  and  slight  evidence  for  a  pre-event  rise  in  Jz 
and  k(Jz).  The  flare-productive  regions  have  both  a  larger  magnitude  of  the  twist  parameter  a 
and  a  greater  variation  over  the  region,  with  a  larger  error  (uncertainty)  in  the  as  parameter 
for  both  flare-productive  regions.  The  magnetic  helicity  density  and  its  parameterizations  show 
a  number  of  unique  flare-event  signatures,  including  significant  Hc>net  and  hc  for  the  immediately 
flare-productive  regions,  and  a  clear  decrease  in  the  latter  prior  to  two  of  the  three  events  detailed 
in  [18].  We  do  find  that  magnetic  shear  was  larger  and  more  widespread  in  the  flare-producing 
regions. 

When  pre-event  variations  and  levels  are  examined  using  both  flare-producing  and  flare-quiet 
regions,  the  obvious,  and  event-unique  signatures  are  rare.  We  do  focus  on  only  three  regions  and 
as  many  flare  events  here;  in  [19]  (see  §4,  below)  we  take  an  in-depth,  statistical  look  at  the  question 
of  flare- unique  signatures. 

4  The  Discriminant  Function  Algorithm  as  applied  to  Flare  Prediction: 

Leka  &  Barnes  Astrophys.  J.,  2003b. 

4.1  Abstract 

Using  Imaging  Vector  Magnetograph  observations  of  photospheric  vector  magnetic  fields  in  seven 
active  regions,  we  apply  discriminant  analysis  to  the  wide  range  of  parameters  discussed  in  Paper 
I  (Leka  &  Barnes  2003,  cf.  §3,  above).  We  treat  the  mean  value  and  the  rate  of  change  of  each 
parameter  as  a  separate  variable.  In  order  to  discover  which  of  these  properties  are  associated 
with  flaring,  we  begin  by  considering  pairs  of  variables.  Applying  Hotelling’s  T2-test  to  the  data, 
we  find  that  a  number  of  properties  show  statistical  differences  between  flaring  and  flare-quiet 
regions.  However,  discriminant  functions  of  pairs  of  variables  have  high  error  rates,  implying  a 
large  degree  of  overlap  in  the  properties  of  flaring  and  flare-quiet  regions.  To  reduce  this  overlap, 
we  simultaneously  consider  larger  numbers  of  variables.  This  results  in  much  lower  estimates  for 
the  error  rate,  but  we  can  find  no  unique  combination  of  variables  which  is  significantly  better 
than  other  combinations.  Our  sample  sizes  are  too  small  for  us  to  directly  compare  the  predictive 
power  of  large  numbers  of  variables,  to  determine  which  are  the  most  important  in  flaring.  Instead, 
we  consider  discriminant  functions  of  all  possible  permutations  of  four  variables.  We  rank  the 
permutations  based  on  the  probability  that  the  samples  come  from  different  populations,  and  look 
for  the  most  frequently  appearing  variables  in  the  best  permutations.  We  present  these  variables 
as  the  most  likely  to  be  associated  with  flaring.  Individually,  these  variables  may  have  little  ability 
to  distinguish  between  flaring  and  flare-quiet  regions,  so  it  is  important  to  consider  combinations 
of  several  variables.  We  construct  a  single  discriminant  function  of  ten  variables  to  support  our 
method  of  selecting  the  most  frequently  occurring  variables  as  being  the  most  important  in  flaring. 
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Table  1:  Parameters  used  in  the  discriminant  analysis 


variable 


8  =  median(AJ) 


vertical  magnetic  field 

B~z 

o{Bz) 

<=(Bz) 

k(Bz) 

=  E\B*\dA 
|«n6t|  =  IE-B*<W| 

horizontal  magnetic  field 

B~h 

a(Bh) 

<Bh) 

KjM _ 


description _ 


seeing,  t.e.,  median  of  the  granulation  contrast 
Distribution  of  Magnetic  Fields  _ 


Bz  =  Bez 

mean  of  vertical  magnetic  field 

standard  deviation  of  vertical  magnetic  field 

skew  of  vertical  magnetic  field 

kurtosis  of  vertical  magnetic  field 

total  unsigned  flux 

absolute  value  of  the  net  flux 

Bh  =  y/B*  +  B  l 

mean  of  horizontal  magnetic  field 

standard  deviation  of  horizontal  magnetic  field 

skew  of  horizontal  magnetic  field 

kurtosis  of  horizontal  magnetic  field _ 


Distribution  of  Magnetic  Field  Spatial  Gradients 


gradient  of  total  field 

\m 

cr(\VhB\) 

<(\VhB\) 

«(|VhB|) 

gradient  of  vertical  field 

WM 

<r(\VhBz\) 

K(\VhBz\) 

gradient  of  horizontal  field 

iv^i 

a(\VhBh\) 

«]VhBh\) 

K(\VhBh\) 


vertical  current  density 

Tz 

a{Jz) 

c(Jz) 

k(Jz) 

Itot  =  Z\Jz\dA 
Inet  =  \EJzdA\ 

1%'t  =  I E  UBZ  >  0)  dA\  +  I  £  JZ(BZ  <  0)  dA\ 
vertical  component  of  hetergeneity  current  density[33] 

JK 

<r(fo 

K(Jz) 

lL  =  ZW?\dA 

It  =  I E  JzhdA\ _ _ 


\VhB\  =  ^{dB/dxf  +  (OB/dy) 
mean  of  gradient  of  field 
standard  deviation  of  gradient  of  field 
skew  of  gradient  of  field 
kurtosis  of  gradient  of  field 

\vhBz\  =  s/WJW  +  W7IW 

mean  of  gradient  of  vertical  field 

standard  deviation  of  gradient  of  vertical  field 

skew  of  gradient  of  vertical  field 

kurtosis  of  gradient  of  vertical  field 

|VfcRh|  =  s/{dBh/dxy  +  (dBh/dy)2 

mean  of  gradient  of  horizontal  field 

standard  deviation  of  gradient  of  horizontal  field 

skew  of  gradient  of  horizontal  field 

kurtosis  of  gradient  of  horizontal  field _ 


Vertical  Current  Density  _ 


Jz  =  ( dBy/dx  -  dBx/dy)/fio 
mean  of  vertical  current  density 
standard  deviation  of  vertical  current  density 
skew  of  vertical  current  density 
kurtosis  of  vertical  current  density 
total  unsigned  current 
absolute  value  of  the  net  current 
[))  dA\  sum  of  absolute  value  of  net  currents  in  each  polarity 
msity[33]  Jz  =  ( bvdBx/dy  —  bxdBy/dx)/fXo 

mean  of  vertical  hetergeneity  current  density 
standard  deviation  of  vertical  hetergeneity  current  density 
skew  of  vertical  hetergeneity  current  density 
kurtosis  of  vertical  hetergeneity  current  density 
total  unsigned  hetergeneity  current 

absolute  value  of  the  net  hetergeneity  current _ 
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Table  1:  -  Continued 


variable 

description 

Current  Helicity 

current  helicity[2] 

hc  =  CBz{dBy/dx  -  dBxjdy ) 

hc 

mean  of  current  helicity 

a(hc) 

standard  deviation  of  current  helicity 

<;(hc) 

skew  of  current  helicity 

K(hc) 

kurtosis  of  current  helicity 

Hc,tot  =  E\hcdA\ 

total  unsigned  current  helicity 

\HCiTiet\  =  | 52  hc  dA\ 

absolute  value  of  net  current  helicity 

Inclination  Angle 

inclination  andle 

7  =  tan  l(Bz/Bh) 

7 

mean  of  inclination  angle 

<7(7) 

standard  deviation  of  inclination  angle 

c(t) 

skew  of  inclination  angle 

«( 7) 

kurtosis  of  inclination  angle 

Shear  Angle 

3-D  Shear  Angle[32] 

=  cos-1  {BP -B0/ BP  B°) 

mean  of  shear  angle 

<*(*) 

standard  deviation  of  shear  angle 

skew  of  shear  angle 

«(*) 

kurtosis  of  shear  angle 

>  45°) 

area  with  shear  >  45° 

A(y  >  80°) 

area  with  shear  >  80° 

Vnl 

mean  of  neutral  line  shear  angle 

cr^NIy) 

standard  deviation  of  neutral  line  shear  angle 

s(Vnl) 

skew  of  neutral  line  shear  angle 

k^nl) 

kurtosis  of  neutral  line  shear  angle 

Wnl  >  45°) 

length  of  neutral  line  with  shear  >45° 

L{Vnl  >  80°) 

length  of  neutral  line  with  shear  >  80° 

Horizontal  Shear  Angle 

^  =  cos  ~l(Bl-B°h/BphB°h) 

mean  of  horizontal  shear  angle 

a(ip) 

standard  deviation  of  horizontal  shear  angle 

skew  of  horizontal  shear  angle 

*#) 

kurtosis  of  horizontal  shear  angle 

A{tl>  >  45°) 

area  with  horizontal  shear  >  45° 

A(V>  >  80°) 

area  with  horizontal  shear  >  80° 

IpNL 

mean  of  neutral  line  horizontal  shear  angle 

standard  deviation  of  neutral  line  horizontal  shear  angle 

silpNL) 

skew  of  neutral  line  horizontal  shear  angle 

kurtosis  of  neutral  line  horizontal  shear  angle 

Wnl  >  45°) 

length  of  neutral  line  with  horizontal  shear  >  45° 

L(i>NL  >  80°) 

length  of  neutral  line  with  horizontal  shear  >  80° 
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Table  1:  -  Continued 


variable 

description 

Twist  Parameter 

twist  parameter 

a  =  CJz/Bz 

a 

mean  of  twist 

a(a) 

standard  deviation  of  twist 

<(«) 

skew  of  twist 

«(«) 

kurtosis  of  twist 

best  fit  aff[21] 

B  =  affVxB 

l«ffl 

absolute  value  of  best  fit  as 

Pseudo- Free  Energy 

pseudo-free  energy  density[32] 

Pf  =  (Bp  —  B°)  z/8ir 

Pf 

mean  of  free  energy 

<r(p/) 

standard  deviation  of  free  energy 

c(P/) 

skew  of  free  energy 

«(p/) 

kurtosis  of  free  energy 

Ef  -  Y,PfdA 

toted  free  energy 

For  each  of  these  parameters,  we  consider  the  mean  value  for  an  epoch,  denoted  by  ()  and  the 
slope  of  a  regression  line,  denoted  by  dj dt. 
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4.2  Discriminant  Analysis 

Our  problem  is  essentially  one  of  statistical  discrimination:  measurements  of  a  number  of  parame¬ 
ters  derived  from  the  photospheric  vector  field  data  (cf.  §3)  come  from  samples  of  two  populations: 
flare  events  or  flare-quiet  events.  The  goal  is  to  determine  whether  the  two  populations  are  statis¬ 
tically  different  from  each  other,  based  on  those  measurements.  If  they  are  statistically  different 
from  each  other,  then  given  a  new  data  point  the  question  is  how  to  classify  it  as  belonging  to 
one  group  of  the  other.  The  differences  are  evaluated  by  performing  Hotelling’s  T2-test,  and  we 
construct  discriminant  functions  in  order  to  predict  group  membership  [11,  1].  The  discriminant 
function  is  constructed  to  maximize  correct  predictions  while  at  the  same  time  maintaining  equal 
probabilities  for  predicting  that  a  region  will  flare  when  it  does  not,  and  for  predicting  that  a  region 
will  not  flare  when  it  does. 

The  most  straightforward  way  to  test  how  well  the  function  discriminates  between  flaring  and 
non-flaring  regions  is  to  construct  a  classification  or  “truth”  table.  Unfortunately,  it  can  be  shown 
that  this  method  always  underestimates  the  number  of  incorrect  predictions  [10].  Other  methods 
can  place  stricter  bounds  on  the  probability  of  misclassification.  For  example,  constructing  a 
discriminant  function  with  one  data  point  left  out,  and  using  that  discriminant  function  to  classify 
the  excluded  point,  then  repeating  this  for  all  data  points  gives  an  unbiased  estimate  of  the  error 
rate  for  the  discriminant  function  with  n  -  1  data  points.  In  the  limit  of  large  n,  this  will  be  an 
estimate  of  the  actual  error  rate  for  the  n  data  point  case.  Our  sample  sizes  are  sufficiently  small 
that  we  restrict  our  use  of  this  approach  to  illustration  only. 

The  limitations  of  these  procedures  include  the  following.  A  high  probability  that  the  samples 
come  from  different  populations  does  not  necessarily  mean  that  the  discriminant  function  will  have 
a  low  error  rate  in  predicting  membership  for  future  observations.  For  example,  two  samples  may 
have  significantly  different  but  also  a  significant  scatter  such  that  there  results  a  great  deal  of 
overlap  in  parameter  space.  A  discriminant  function  will  never  have  a  low  error  rate  for  such  a 
situation.  However,  if  the  actual  error  rate  can  be  accurately  estimated,  then  the  combination  of 
the  significance  of  the  difference  in  the  samples  combined  with  the  estimate  of  the  error  rate  can 
be  used  to  learn  something  about  the  samples.  If  the  samples  are  significantly  different  and  the 
actual  error  rate  is  low,  then  the  samples  occupy  distinct  regions  in  the  space  of  the  variables  used 
in  the  discriminant  function.  This  is  the  ideal  case,  in  which  it  is  possible  to  definitively  describe 
those  physical  properties  of  an  active  region  sufficient  and  necessary  for  it  to  flare.  If  the  error  rate 
is  high,  despite  a  significant  difference  between  the  samples,  then  there  is  a  great  deal  of  overlap 


Table  2:  Active  Regions  used  in  Discriminant  Analysis 


NO  A  A  AR 
number 

Date 

Coordinates 

Area 

(^H) 

Mag. 

class 

McIntosh 

class 

Observed  GOES 
SXR  flares 

8210 

1998  May  01 

S18  W05 

270 

07 

Cho 

C2.8,  C2.6,  M1.2 

8636 

1999  Jul  23 

N20  E04 

550 

07  6 

Fki 

Ml.l 

8771 

1999  Nov  25 

S15  W48 

750 

07 

Eki 

C1.6,  M2.0 

8891 

2000  Mar  01 

S15  Ell 

1030 

07 

Eki 

•  .  . 

9026 

2000  Jun  05 

N22  E20 

870 

07<* 

Fki 

C3.8 

9165 

2000  Sep  15 

N13  E03 

140 

07 

Dai 

C7.4 

0030 

2002  Jul  15 

N18  E03 

781 

07  8 

Fkc 

X3.0,  M1.8 
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Figure  5:  The  discriminant  function  for  the  variables  dn(Jz)/dt].  Flaring  epochs  (+)  are  color 

coded  for  the  Soft  X-Ray  flare  class,  C  (green)  M  (orange)  and  X  (red);  quiet  epochs  are  also  indicated  (o). 
Each  epoch  is  labeled  with  its  parent  NOAA  AR  number.  The  means  of  each  sample  are  shown  with  circles, 
and  the  resulting  discriminant  function  is  the  solid  line. 


in  the  parameter  space.  This  suggests  that  there  is  an  additional  physical  property  necessary  to 
produce  an  energetic  event.  In  order  to  make  either  of  the  above  statements,  an  accurate  estimate 
of  the  actual  error  rate  must  be  available  which  is  not  generally  the  case.  The  number  of  incorrect 
predictions  in  the  classification  table  provides  a  lower  bound,  and  better  estimates  can  be  made, 
particularly  for  large  sample  sizes,  but  the  actual  error  rate  will  not  be  known. 

4.3  Application  to  IVM  Data 

For  this  study,  Imaging  Vector  Magnetograph  data  for  all  seven  active  regions  (see  Table  2)  are 
included.  Each  time  sequence  of  magnetograms  is  divided  into  different  epochs,  based  on  the 
occurrences  of  flares  and  on  tune  gaps  in  the  data.  For  this  demonstration,  we  have  a  group  of  10 
flaring  epochs,  and  a  group  of  14  non-flaring  epochs. 

The  behavior  of  each  parameter  in  an  epoch  is  described  by  its  weighted  mean,  i.e.,  the  param¬ 
eter’s  overall  magnitude,  and  the  slope  of  a  weighted  regression  line  fit  to  the  temporally  sampled 
data,  i.e.,  the  temporal  rate  of  change  of  the  parameter  in  question.  A  complete  list  of  the  param¬ 
eters  considered  is  given  in  Table  1. 

4.3.1  A  Two  Variable  Discriminant  Function 

We  first  consider  the  combination  of  the  mean  of  the  skew  of  the  current  helicity,  fc(hc))  and 
the  slope  of  the  kurtosis  of  the  vertical  current,  dK,(Jz)/dt].  Using  these  two  variables  alone,  the 
probability  that  the  samples  are  from  distinct  populations  is  0.95,  the  highest  probability  of  any 
variable  pair.  The  discriminant  function  for  this  case  is 

=  0.12  -1.10fc(M>  +  L03^jgp  (»> 

in  standardized  variables.  The  nearly  equal  magnitudes  of  the  resulting  coefl&cients  implies  that 
they  have  approximately  equal  predictive  power.  This  function,  along  with  the  data  from  which  it 
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was  constructed,  is  shown  in  Figure  5.  Note  that  in  this  case,  there  is  neither  obvious  clustering  of 
different  epochs  from  the  same  active  region,  nor  of  the  flare  classes. 


Table  3:  2- Variable  Discriminant  Function  Classification  Table 

predicted 


observed 


flare 

no  flare 

flare 

7 

3 

no  flare 

4 

10 

The  classification  table  for  this  discriminant  function  is  shown  in  Table  3;  from  either  the 
table,  or  by  directly  inspecting  the  figure,  one  finds  that  seven  of  twenty-four  classifications  were 
made  incorrectly,  for  an  error  rate  of  0.29.  The  “n  -  1  data  points”  method  (outlined  above) 
results  in  an  estimated  error  rate  identical  to  the  table.  Hence,  despite  the  high  probability  for 
distinct  populations,  the  high  estimates  for  the  error  rate  suggest  that  these  two  variables  alone 
are  insufficient  to  clearly  distinguish  flaring  from  flare-quiet  epochs. 

4.3.2  Six  Variable  Discriminant  Functions 

To  improve  upon  the  error  rate  for  predicting  the  membership  of  a  new  observation,  the  analysis 
must  be  performed  by  considering  more  variables  simultaneously.  Ideally,  a  single  discriminant 
function  is  constructed  of  all  the  variables,  and  the  relative  magnitudes  of  the  coefficients  are  used 
to  estimate  the  predictive  power  of  each  parameter  relative  to  the  others.  However,  due  to  our 
small  sample  sizes  we  are  currently  unable  to  consider  more  than  a  few  variables  simultaneously. 
Still,  by  considering  all  the  possible  permutations  of  parameters  for  a  given  number  of  variables, 
we  can  sort  the  statistical  results  to  determine  the  best  discriminant  function  of  that  number  of 
variables. 

Doing  this  for  the  four-variable  case,  we  take  the  variables  appearing  most  frequently  in  the 
best  discriminant  functions  to  be  the  parameters  with  the  most  predictive  power  (this  approach  for 
greater  than  four  variables  becomes  prohibitive  since  the  number  of  permutations  grows  so  rapidly). 
Nevertheless,  using  the  best  four-variable  cases  as  starting  points,  we  are  able  to  arrive  at  several 
interesting  examples  of  discriminant  functions  of  six  variables.  One  such  example  is 

/  =  1.2  +  31.0(|V^|>  +  6.9<c(|  VfcB&|))  +  10.4</,fet>  +  81.5<*)  +  66.6(?M)  +  7.2  ^"^(4) 

in  standardized  variables  (see  Table  1  for  variable  descriptions).  The  probability  that  flaring  and 
flare-quiet  epochs  come  from  different  populations  for  these  variables  is  0.999999,  and  the  discrim¬ 
inant  function  is  able  to  successfully  classify  all  of  the  data  points,  i.e.,  the  error  rate  is  0.0  (see 
Table  4).  However,  this  is  an  example  of  using  a  truth  table  with  caution:  the  “n  —  1  data  points” 
error-estimation  method  gives  an  error  rate  estimated  at  0.083. 

A  second  interesting  example  the  discriminant  function 

/  =  1.0  -  38.8 ($ tot)  +  6.8(<(|  VhBh\)}  +  36.8(Jtot)  +  27.6(<r(tf ))  +  33.6(k(^))  +  2.1rf|^wet|(5) 

in  standardized  variables,  which  results  in  a  different-population  probability  of  0.999996,  and  the 
discriminant  function  mis-classifies  one  data  point  for  an  error  rate  of  0.042  (see  Table  4),  identical 
to  the  estimated  actual  error  rate  derived  from  the  “n  —  1”  approach. 
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Table  4:  6- Variable  Discriminant  Function  Classification  Tables 

predicted  predicted 


observed 


flare 

no  flare 

flare 

10 

0 

no  flare 

0 

14 

observed 


flare 

no  flare 

flare 

9 

1 

no  flare 

0 

14 

Evidently,  there  exist  combinations  of  variables  for  which  flaring  and  flare-quiet  epochs  occupy 
quite  different  regions  in  phase  space.  The  fact  that  there  are  only  two  variables  which  appear 
in  both  examples  above  suggests  that  discriminant  functions  of  larger  numbers  of  variables  would 
be  preferable,  since  there  are  distinctions  between  the  groups  which  are  not  being  made  in  each 
example.  However,  to  do  this  simply  requires  more  data. 

4.4  Summary 

While  flaring  and  flare-quiet  epochs  appear  to  be  statistically  different  for  groups  of  variables 
or  even  some  pairs  of  variables,  at  this  time  we  find  no  single,  or  even  small  number  of,  physical 
properties  of  an  active  region  which  directly  results  in  flaring.  We  do  demonstrate,  however,  that  to 
accurately  predict  which  active  regions  will  flare,  many  variables  must  be  considered  simultaneously. 
By  considering  all  permutations  of  four  variable  discriminant  functions,  we  conclude  that  the  active 
regions  which  flare  are  likely  to  have:  an  increasing  kurtosis  of  the  vertical  current,  dn(Jz)/dt,  an 
increasing  kurtosis  of  the  twist  parameter,  dn(a)/dt,  a  decreasing  kurtosis  of  the  current  helicity, 
dn{hc)/dt ,  and  a  large  mean  and  standard  deviation  of  the  neutral  line  horizontal  shear  angle, 
$nl)  and  (<t{$nl))-  The  kurtosis  of  a  distribution  describes  the  non-normal  appearance  of  either 
of  a  distribution’s  wings;  thus  these  properties  roughly  correspond  to  having  an  increase  in  the 
magnitudes  of  the  strongest  vertical  currents,  an  increase  in  the  strongest  twist  parameters,  a 
decrease  in  the  strongest  magnitudes  of  the  current  helicity  and  a  neutral  line  along  which  the  field 
is  highly  non-potential.  Other  properties  are  likely  to  be  necessary  for  flaring,  but  these  are  the 
most  frequently  occurring  properties  in  our  best  four  variable  discriminant  functions. 

Our  six  variable  examples  indicate  that  while  combinations  can  be  found  which  result  in  quite 
good  predictions,  the  combinations  are  not  unique  and  hence  larger  numbers  of  variables  need  to  be 
considered  simultaneously.  To  do  this  requires  more  data  in  order  to  simply  consider  more  variables 
and  also  to  make  accurate  estimates  of  the  error  rates. 

5  Quantifying  the  Coronal  Topology  with  the  Magnetic  Charge  Topology 

Barnes,  Leka  &:  Longcope  Astrophys.  J.,  2003. 

One  of  the  great  difficulties  encountered  in  trying  to  understand  solar  flares  is  that  one  must  gener¬ 
ally  understand  the  magnetic  field  in  the  corona,  where  the  flare-onset  energy  release  occurs,  using 
extrapolations  from  measurements  of  vector  magnetic  fields  at  the  photosphere.  The  Minimum 
Current  Corona  (MCC)  model  [24]  is  one  approach  which  has  shown  some  success.  This  model 
has  been  applied  to  a  single  magnetogram  of  NOAA  AR  6993/6994  [26],  and  was  able  to  predict 
some  of  the  field  lines  along  which  the  flares  appeared  to  occur.  The  MCC  is  an  extension  of  the 
Magnetic  Charge  Topology  (MCT)  models  discussed  by  numerous  authors  [3,  8,  30,  16,  6,  29]. 

In  the  MCT  models,  the  photospheric  Bz  field  is  used  to  derive  the  locations  and  strengths 
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-150  -100  -50  0  50  -150  -100  -50  0  50 

Figure  6:  One  magnetogram  from  NOAA  AR  8636,  showing  the  smoothed  field,  the  pole  locations,  the 
boundaries  of  the  partitions  defining  the  pole  locations  (black  lines)  and  the  connections  between  poles  (blue 
lines)  in  the  left  panel,  and  the  locations  of  the  poles  and  nulls  (blue,  red  and  orange),  the  null  spines  (solid 
black  lines),  the  fan  traces  (dotted  black  lines)  and  some  of  the  separators  (green  lines)  in  the  right  panel. 


of  a  set  of  magnetic  monopoles;  the  field  everywhere  above  the  photosphere  is  then  assumed  to 
be  due  only  to  these  monopoles.  Knowing  the  magnetic  field,  we  can  in  principle  determine  the 
locations  of  magnetic  null  points,  where  the  magnetic  field  vector  vanishes.  In  this  approximation, 
each  field  line  begins  and  ends  on  a  monopole  or  null  (or  at  infinity),  which  provides  a  way  to 
quantify  the  topology  of  the  coronal  magnetic  field.  We  define  the  connectivity  matrix,  V’ij,  as  the 
flux  connecting  the  ith  pole  to  the  jth.  pole.  The  surfaces  which  bound  the  flux  bundles  joining  two 
poles  are  known  as  separatrix  surfaces.  The  intersection  of  two  separatrix  surfaces  is  a  separator, 
which  is  a  field  line  connecting  two  null  points.  It  has  been  shown  [9,  17]  that  separators  are  the 
site  of  magnetic  reconnection,  and  hence  likely  to  be  the  site  of  energy  release  in  flaring. 

5.1  Application  of  MCT  to  IVM  data 

We  have  adapted  the  code  developed  by  Dr.  D.W.  Longcope  to  apply  the  MCT  model  to  time- 
series  of  IVM  data  for  several  active  regions  (cf.  Table  2).  The  code  had  previously  been  applied 
to  a  single  lower  spatial-resolution  magnetogram  of  a  simple  active  region;  as  such,  it  has  required 
extensive  modification  to  be  used  with  our  datasets.  First,  the  magnetogram  is  smoothed,  using 
a  potential  field  extrapolation;  only  areas  of  the  magnetograms  in  which  the  measured  vertical 
field  strength  is  not  consistent  with  zero  are  used,  determined  by  the  measured  noise  levels.  The 
remaining  part  of  the  magnetogram  is  partitioned  using  a  downhill  gradient  method  where  every 
pixel  is  associated  with  its  local  maximum  or  minimum.  This  results  in  such  a  large  number  of 
partitions  as  to  be  unwieldy,  so  partitions  with  small  flux  are  merged  with  their  closest  neighbor 
of  the  same  polarity,  provided  that  neighbor  is  within  a  maximum  distance.  Our  motivation  in 
doing  this  is  to  simplify  areas  of  plage,  where  small  variations  in  field  strength  over  large  areas  can 
result  in  many  small  partitions,  while  retaining  strong  partitions  with  polarity  opposite  to  their 
surroundings,  such  as  delta  spots.  Each  final  partition  is  then  represented  by  a  magnetic  pole,  the 
strength  of  which  is  determined  by  the  total  flux  of  the  partition,  and  the  location  of  which  is  the 
flux- weighted  average  position  of  the  partition  (Figure  6). 

The  connectivity  matrix  is  calculated  using  a  Monte-Carlo  method  to  select  field  lines  emanating 
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log  flux  in  a  connection 


Figure  7:  Panel  (a)  shows  the  number  of  connections  found  as  a  function  of  the  number  of  field  lines  used 
in  the  calculation.  Panel  (b)  plots  the  number  of  connections  found  as  a  function  of  the  log  of  the  flux  in 
the  connections  for  the  case  of  2  x  10s  field  lines  (solid  line)  and  the  case  of  2  x  104  field  lines  (dotted  line). 


from  one  pole  and  tracing  them  to  their  termination  on  another  pole  (or  at  infinity).  In  the  initial 
application  of  the  code,  300  field  lines  originating  on  the  positive  poles,  and  an  equal  number 
originating  on  the  negative  poles,  were  sufficient  to  calculate  the  connectivity.  In  applying  it  to 
IVM  data,  we  found  300  to  be  severely  inadequate  (Figure  7).  We  have  modified  Dr.  Longcope’s 
initial  approach  to  use  at  least  2  x  104  field  lines,  with  the  fraction  of  these  field  lines  originating 
from  a  given  pole  determined  by  the  fraction  of  the  total  flux  it  contains,  with  a  minimum  of  10 
field  lines.  In  Figure  7(a)  we  show  the  number  of  connections  found  as  a  function  of  the  number  of 
field  lines  used  in  the  calculation  for  one  magnetogram;  N.B.,  the  result  for  300  is  not  even  shown. 
While  we  risk  missing  some  connections  even  with  2  x  104  field  lines,  using  a  significantly  higher 
number  is  prohibitively  time  consuming.  Additionally,  we  have  determined  that  the  connections 
which  are  missed  are  preferentially  those  with  very  small  fluxes  [see  Figure  7(b)],  which  are  unlikely 
to  contribute  significantly  to  flaring  activity. 

In  addition  to  the  connectivity  matrix,  we  calculate  the  locations  of  null  points  and  separa¬ 
tors.  The  original  null-finding  procedure  from  Dr.  Longcope  was  unable  to  find  all  of  the  nulls 
in  complicated  pole  arrangements  such  as  those  we  presented.  We  have  improved  upon  the  algo¬ 
rithm,  but  in  essence  the  problem  is  one  of  multi-dimensional  root  finding,  for  which  there  exists 
no  general  solution.  Likewise,  tracing  the  separators  depends  on  following  a  single  field  line  from 
which  neighboring  field  lines  are  exponentially  diverging.  In  some  cases,  this  can  be  accomplished, 
but  not  in  general.  Figure  6  shows,  as  a  figure,  the  connectivity  matrix  ipij,  the  locations  of  the 
nulls  and  of  some  separators.  Note  that  all  these  quantities  are  not  independent.  When  only  prone 
photospheric  nulls  are  present  [25], 

Nsep  =  Nc  —  Nn  —  1  (6) 

where  Nsep  is  the  number  of  separators  above  the  photosphere,  Nc  is  the  number  of  connections 
and  Nn  is  the  number  of  nulls.  Our  MCT  models  of  IVM  data  contain  other  types  of  nulls  (upright 
nulls),  which  can  also  result  in  separators  lying  in  the  plane  of  the  photosphere.  In  conjunction 
with  Dr.  Longcope,  we  hope  to  expand  this  treatment  to  allow  for  such  nulls  and  separators  to 
more  adequately  represent  the  observations. 
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6.  Summary  of  Accomplishments 


19.0  19.5  20.0  20.5  21.0 
UTtime 


17.0  17.5  18.0  18.5  19.0 
UTtime 


[  *********  ) 

[  8891  J 

20.0  20.5  21.0  21.5 
UTtime 


Figure  8:  The  time  evolution  of  the  number  of  separators  for  NOAA  AR  0030, 8636  and  8891.  The  red  line 
indicates  the  time  of  an  X-class  flare,  while  the  yellow  lines  indicate  the  times  of  M-class  flares. 


5.2  Quantifying  the  Coronal  Topology:  Sample  Results 

In  essence,  this  is  the  foundation  for  the  MCT  model;  evolving  this  model  in  time  leads  to  currents 
along  separators,  which  is  the  origin  of  the  term  Minimum  Current  Corona.  However,  even  before 
considering  the  currents,  we  can  quantify  the  complexity  of  the  coronal  magnetic  field  using  the 
information  contained  in  the  connectivity  matrix  and  related  quantities.  For  example,  since  recon¬ 
nection  is  likely  to  occur  at  the  intersection  of  separatrix  surfaces,  we  consider  the  evolution  of  the 
number  of  separators  (Figure  8),  as  determined  by  equation  (6).  Three  active  regions  are  studied 
in  detail,  two  of  which  flared  during  the  period  of  observation  (NOAA  AR  8636,  0030),  while  the 
third  did  not  (NOAA  AR  8891).  Although  there  is  significant  scatter  in  the  number  of  separators, 
the  regions  which  flared  do  generally  have  significantly  more  separators  than  the  region  which  did 
not  flare. 

The  MCT  model  cam  be  used  to  quantitatively  characterize  the  properties  of  the  coronal  mag¬ 
netic  field,  rather  than  the  directly  measured  photospheric  field.  In  the  context  of  flare  prediction, 
we  can  consider  three  broad  categories  of  possible  flare-predicting  parameters,  with  some  overlap. 
The  first  category  simply  measures  the  overall  complexity  of  a  region.  If  the  reconnection  which 
powers  a  flare  is  a  means  of  simplifying  the  topology  of  the  coronal  field,  then  presumably  the 
field  must  start  with  a  certain  amount  of  complexity,  and  more  complex  regions  are  more  likely 
to  flare.  The  second  category  looks  for  a  trigger  for  a  flare.  For  example,  the  emergence  of  new 
flux  into  an  existing  active  regions  is  thought  to  be  a  mechanism  for  triggering  a  flare.  Finally, 
if  flares  are  indeed  associated  with  reconnection  along  a  separator,  then  such  separators  must  be 
present  in  an  active  region,  and  we  can  consider  the  number  of  possible  locations  at  which  a  flare 
fAn  occur.  Whether  parameters  describing  any  of  these  categories  will  prove  useful  in  flare  predic¬ 
tion  must  be  tested  on  a  statistical  basis.  With  the  discriminant  analysis  already  applied  to  the 
photospheric  parameters  (cf.  §4),  we  have  the  machinery  in  place  to  make  such  a  test  of  inferred 
coronal  parameters. 

6  Summary  of  Accomplishments 

To  summarize,  in  the  last  three  years  of  effort  under  this  contract,  we  have: 

•  Acquired,  reduced  and  prepared  IVM  data  for  six  active  regions  for  analysis;  this  includes 
selecting  target  regions  that  included  flare-productive  and,  as  counter-examples,  flare-quiet 
epochs. 
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6.  Summary  of  Accomplishments 


•  Developed  a  successful  method  for  access,  storage,  and  manipulation  of  large  magnetogram 
data-sets  using  a  platform-independent  approach. 

•  Developed  an  algorithm  to  model,  and  account  for,  the  seeing-induced  variations  in  the 
observed  vector  magnetic  field  data  and  quantities  derived  from  it. 

•  Derived  and  examined  the  magnitude  and  variation  of  more  than  80  parameters  that  describe 
the  magnetic  state  of  the  solar  photosphere  as  it  relates  to  flare  events. 

•  Developed  a  statistical  approach  based  on  Discriminant  Function  Analysis  and  Hotelling’s 
T2-test,  to  evaluate  the  possibility  of  distinguishing  between  flare-imminent  and  flare-quiet 
active  regions  and  epochs  using  the  parameterizations  derived  from  time-sequences  of  vector 
magnetic  field  data. 

•  Applied  the  statistical  tests  to  the  24  flare-productive  and  flare-quiet  epochs  available,  de¬ 
riving  preliminary  lists  of  parameters  which  may  be  useful  to  both  the  flare-prediction  and  the 
modeling  communities. 

•  Attempted  to  apply  the  Magnetic  Charge  Topology  method  to  quantifying  the  coronal  topol¬ 
ogy  and  its  complexity.  Upon  discovery  that  the  codes  as  delivered  to  us  from  Dr.  Longcope 
were  not  suited  for  data  as  complex  as  the  flaring  active  regions  require,  we  have  embarked 
upon  substantial  modifications  to  the  codes  (which  are  being  developed  in  collaboration  with 
Dr.  Longcope  and  shared  with  him). 

•  Nonetheless,  we  have  presented  a  preliminary  example  of  the  temporal  evolution  of  the 
coronal  complexity  in  three  active  regions. 
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ATTACHMENT 2:  CODELIST 


Table  1:  IVM  Data  Handling  and  Parametrization  IDL  Code  List1 


analyze.jjstruc.pro 

calc_errors.pro 

cross  .correl.  pro 

define.jjstruc.pro 

do.divBh.pro 

do_run.plot_params.pro 

expand.point.pro 

getimg_mm.pro 

ivm_mdi_align.pro 

loadofFtrim.pro 

moment.err.pro 

outofres.pro 

plot_flare_params.pro 

read.jjstruc.pro 

run_plot.params.pro 

see.data.pro 

startit.pro 

timetostring.pro 

yrdaytodate.pro _ 


apply.offset.pro 

ch.pointinfo.pro 

csh.jj.pro 

derive_seeing.pro 

doJix_analyze.jjstruc.pro 

doublebz.pro 

extract  .date,  pro 

get.offset.pro 

jjivm.mask.pro 

make.param.pro 

moment  .weight .  pro 

paper.plot  Jcimage.pro 

plot.qNd.pro 

read.offsets.pro 

save.jjtime.pro 

select_data.box.pro 

test.bperr.pro 

trim_array.pro 

zero.entries.pro 


avg.times.pro 

converttime.pro 

cshear.jjstruc.pro 

do_analyze.jjstruc.pro 

do_paper_plot.params.pro 

doy.pro 

fixshear.pro 

intores.pro 

load_analyze.pro 

mask.pukas.pro 

my_shift.pro 

paper  _run_plot  .par  ams .  pro 

read.jj_data.pro 

redoquicklooks .  pro 

save.offsets.pro 

shift.center.pro 

time_startstop.pro 

trim_resize.pro 


1  The  IVM  Data  Reduction  code  is  available  at: 

www .  solar .  if  a .  hawaii .  edu/Ref  erence/IVM/IVMjdatajred .  html 


Table  2:  Discriminant  Function  IDL  Code  List 


all.dist.pro 

analyze.dfa.pro 

count.dfa.pro 

dfa.pro 

dfaJake.pro 

dfaJortran.pro 

dist.dfa.pro 

err.write.pro 

inc.pro 

mkflarect.pro 

sortp2.dfa.pro 

sort.dfa.pro 

sort.rd.pro 

Table  3:  Discriminant  Function  FORTRAN  Code  List 


programs 
one_dfa.f 
err.dfa.f 
sort.dfa.f 


subroutines _ 

dfa(ivar,xl,ifll,err) 

ludcmp(a,indx,d) 

lubksb(a,indx,b) 

rd.dfa 


functions 

betacf(a,b,x) 

betai(a,b,x) 

gammln(xx) 


Table  4:  MCT  IDL  Code  List 


addmask.pro 

ad  vance.difLrot .  pro 

all_nulls.pro 

all_seprs.pro 

allmask,  pro 

analyze_mct.pro 

animate-mask.pro 

animate.partition.pro 

ar_err.struc.pro 

arr2mg.pro 

assemble. pro 

assemblel.pro 

assemble-.err.pro 

auto-con.pro 

auto_mask.pro 

auto-maskf.pro 

automate. pro 

automatel.pro 

automatens.pro 

avg_bp.pro 

bipole,  pro 

bp.eval.pro 

calc.err.pro 

clip-fi.pro 

closest-pole.pro 

colorscale.pro 

colorscale2.pro 

compare.pro 

complete-elliptic  Jnt.pro 

connectivity.pro 

correct-con.pro 

correct_conl.pro 

create_null.pro 

create_nu!10.pro 

cribbon.pro 

crit.stress.pro 

ctest.pro 

def_region.pro 

disk2tan_plane.pro 

display.pro 

do-paper_plot-params.pro 

do_report_plot_params.pro 

eigendwl.pro 

ell_func_rc.pro 

ell_func_rd.pro 

ell-func-rf.pro 

err_std.pro 

err-StdO.pro 

err-std.struc.pro 

err-Struc.pro 

evaLa.pro 

eval_a_ff.pro 

eval.b.pro 

eval_b2p.pro 

eval.bJF.pro 

eval_bprime.pro 

eval-bprime_ff.pro 

evaLchi.pro 

evolve_poles.pro 

exp_grid.pro 

filter.pro 

find_all.pro 

find_chrg.pro 

find_chrgl.pro 

find-chrga.pro 

find_null.pro 

fl_from_point.pro 

fl_from_point0.pro 

fljntgrt.pro 

fLview_xform.pro 

fwrite.pns.pro 

gammaJine.pro  * 

get.berrs.pro 

gradient.mask.pro 

heLchar.pro 

imgexp.pro 

imgscUen.pro 

in.region.pro 

init_sepxJine.pro 

interp_poles.pro 

inv  _lam -func.  pr  o 

jjtime2pls.pro 

kpno_rd.pro 

label.pro 

length.vector.pro 

line_info.pro 

line_piece.pro  4 

line-pieceO.pro 

locate_sepr.pro 

make.mask.pro 

mask-file,  pro 

mct-struc.pro 

mdi-rd.pro 

merge_nls.pro 

mergemlsJabel.pro 

merge.pls.pro 

merge.plsl.pro 

merge.pns.pro 

mg.pro 

mgplot.pro 

mkbwt.pro 

mkcolort.pro 

mkflarect.pro 

moment.err.pro 

moment  -weight .  pro 

morejnulls.pro 

most-nulls.pro 

mp.pro 

munpack.pro 

mut  Jnduct  _mat  .pro 

normalize_3dvec.pro 

null  .scan,  pro 

num_pns.pro 

old_run_plot.params.pro 

one.err.pro 

one.err.struc.pro 

paper  _run-plot-params. pro 

par.err.pro 

parameter.pro 

pb_ang.pro 

pixel2tan_plane.pro 

plot -flare-params. pro 

pls-struc.pro 

pns_err_struc.pro 

pns_err_strucl  .pro 

pns_struc.pro 

pns-test.pro 

pole.crash.pro 

post-rx-fl.pro 

randomJieldJine.pro 

rd.connect.pro 

rd_connects.pro 

rd-multi.pro 

rd.nulls.pro 

rd-pns.pro 

rd.poles.pro 

rdnls.pro 

rdpls.pro 

rdsepr.pro 

rd  view. pro 

regchar.pro 

region_anl.pro 

resolve. pro 

run_plot_params.pro 

Table  4:  -  Continued 


scale  _pls.  pro 

seprJine.pro 

sepx_pole_map.pro 

show.3dpoles.pro 

show.alls.pro 

show.con.dev.pro 

show.dev.pro 

showJlux.dist.pro 

show.jjtime.pro 

showmulls.pro 

show.plane.pro 

show.pns_strucl  .pro 

show.sepr.pro 

show_spx.xangles.pro 

spr_polar2xyz.pro 

tan.plane2disk.pro 

track.multi_null.pro 

update  Jiulls.pro 

view_xform.pro 

write.err.pro 

xmgp.pro 

yet_more_nulls.pro 


script.pro 

sepr.plane.pro 

set.brmax.pro 

show_all.pro 

show  .alls.err. pro 

show.connectivity.pro 

show.domain.pro 

show  _flux.tube.pro 

show  Joop  _skelet  on .  pro 

show.one.pmap  .pro 

show.pns.pro 

show.poles.pro 

show  _sepr2.  pro 

show_zO-fl.pro 

spx.xangles.pro 

tan_plane_matrix.pro 

track.null.pro 

upright.pro 

write.connect.pro 

write_pls.pro 

xmgpart.pro 


sepr_info.pro 

sepxJine.pro 

show_3dnulls.pro 

show.all-Struc.pro 

show_ass.pro 

show.cr_skel.pro 

show  .err. pro 

show-gamma  Jines  .pro 

show.mask.pro 

show.partition.pro 

show_pns.struc.pro 

show  _random_3dlines.  pro 

show  jsepx. pro 

solar  _rotate.view.pro 

sun.globe.pro 

testJmg.pro 

triad. pro 

vertical_filter.pro 

write.connects.pro 

write.pns.pro 

xselect_nulls.pro 


Table  5:  MCT  FORTRAN  Code  List 


programs 

subroutines 

functions 

connectivity.f 

eulerstep 

urandJnit 

intflag 

ilev 

cang 

urand 

ran2 

3 


